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Abstract 

The dependence on Q 2 (the negative square of the 4-momentum of the exchanged virtual photon) 
of the generalised Gerasimov-Drell-Hearn integral for the proton has been measured in the range 
1.2 GeV 2 < Q 2 < 12 GeV 2 by scattering longitudinally polarised positrons on a longitudinally 
polarised hydrogen gas target. The contributions of the nucleon-resonance and deep-inelastic regions 
to this integral have been evaluated separately. The latter has been found to dominate for Q 2 > 
3 GeV 2 , while both contributions are important at low Q 2 . The total integral shows no significant 
deviation from a 1 /Q 2 behaviour in the measured Q 2 range, and thus no sign of large effects due to 
either nucleon-resonance excitations or non-leading twist. 
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The Gerasimov-Drell-Hearn (GDH) sum rule Q 
relates the anomalous contribution k in the nu- 
cleon magnetic moment to an energy- weighted in- 
tegral of the difference of the nucleon's total spin- 
dependent photoabsorption cross sections: 

f°° dv 2n 2 a 2 

J vo -*&)}- = - -j^r K - (!) 

where cri^j is the photoabsorption cross section 
for total helicity of the photon-nucleon system equal 
to \ (|), v is the photon energy in the target rest 
frame, vq is the pion production threshold and M 
is the nucleon mass. For the proton (n p = +1.79) 
the GDH sum rule prediction is — 204/xb. The 
importance of this sum rule is due to the fact 
that it is based mostly on very general princi- 
ples of causality, unitarity, crossing symmetry, and 
Lorentz and gauge invariance. It has never been 
directly tested, due to the need for a circularly po- 
larised beam with a longitudinally polarised tar- 
get, and a wide range of photon energies that has 
to be covered. There are available several pre- 
dictions for the contribution of nucleon resonance 
excitation to the GDH integral Q, derived from 
multipole analyses of data for unpolarised single- 
pion photoproduction, and a prediction for the 
contribution of high energy multi-hadron produc- 
tion ||, based on a multiple- Reggeon exchange 
analysis of deep-inelastic asymmetries. The es- 
timate from multipole analysis was confirmed by 
preliminary results from the GDH experiment at 
Mainz Q , which covered the photon energy range 
from 200 MeV up to 800 MeV. 

The GDH integral can be generalised to the 
case of absorption of polarised transverse virtual 
photons with squared four- momentum — Q 2 ff: 

/>oo J 

Igdh(Q 2 )= [ai^Q^-a^Q 2 )}— (2) 

Jo Q 2 A /l + 7 5 

M Jo K x' { ) 

where g\ and g 2 are the polarised structure func- 
tions of the nucleon, 7 2 = Q 2 /v 2 = (2Mx) 2 /Q 2 , 



x = Q 2 /2Mu and x = Q 2 /2Mv . The quantity 
A\ is the longitudinal asymmetry for virtual pho- 
toabsorption, while F\ is the unpolarised structure 
function of the nucleon. The Gilman nota tion || 
for the virtual photon flux factor K = v^J\ + j 2 
has been used. It should be noted that elastic 
scattering occurring at x = 1 does not contribute 
to the generalised integral. Other generalisations 
of the GDH integral also have been considered || . 
They differ from the definition given in Eq. || by 
terms in the integral that are proportional to 7 2 
and which therefore vanish in both the real-photon 
(Q 2 = 0) and the deep-inelastic (Q 2 > 1 GeV 2 
and 7 2 — > 0) limits. Since 7 2 is not small in the 
nucleon-resonance region and at moderate Q 2 (e.g. 
7 2 is larger than unity for the i-33(1232)-resonancc 
for 0.2 GeV 2 < Q 2 < 2 GeV 2 ), these generalisa- 
tions are equivalent for finite values of Q 2 only if 
the contributions of the nucleon-resonance excita- 
tions are small. 

The generalisation of the GDH integral to non- 
zero photon virtuality Q 2 provides a way to study 
the transition from polarised lepton scattering from 
the nucleon, which is dominated by deep-inelastic 
scattering (DIS) at large photon-nucleon centre 
of mass energy W = y 1 M 2 + 2Mv — Q 2 , to the 
polarised real photon absorption on the nucleon, 
which is dominated by nucleon-resonance excita- 
tion at low W. In leading twist (e.g. for Q 2 — * oo), 
Eq. |3| simplifies: the elastic contribution excluded 
from the integral is of higher twist, the factor 
1/yT + 7 2 is 1+C(1/Q 2 ), and in leading twist 
approximation it is generally believed (though not 
formally proven) that the Burkhardt-Cottingham 
sum rule J Q gi{x, Q 2 )dx = holds. This leads to: 

T f^2\ 167r 2 a 

where Ti = J Q g\(x)dx is predicted to have only 
the weak Q 2 -dependence due to QCD evolution. 
As Q 2 — > 0, Igdh(Q 2 ) must change sign in order 
to reach the negative value predicted by the GDH 
sum rule at Q 2 = 0. 

Several phenomenological models have been pro- 
posed to describe the dependence of the gener- 
alised GDH integral on Q 2 f^V^ ■ Some of these 
models predict large effects from nucleon-reson- 
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ance excitation ]|,[T(| or from higher twist pj,|TT|] : 
even for Q 2 up to a few GeV 2 . Other models based 
on chiral perturbation theory have been proposed 



but their application is limited to Q 2 <C 1 GeV 2 12 
The contribution of the region W 2 > 3.24 GeV 2 ^ 
to the GDH integral defined in Eq. was recently 
measured p"3[ ] for the proton and the neutron in 
the range 0.8 GeV 2 < Q 2 < 12 GeV 2 , showing 
that higher-twist effects do not appear to be signif- 
icant in the measured region. This paper presents 
a measurement of the contribution of the reso- 
nance region to the GDH integral for the proton in 
a similar Q 2 -range (1.2 GeV 2 < Q 2 < 12 GeV 2 ). 
In combination with the analysis at higher W 2 , 
this provides the first experimental determination 
of essentially the complete GDH integral for the 
proton over a range of Q 2 values. 

The measurement was performed in 1997 with 
a 27.56 GeV beam of longitudinally polarised po- 
sitrons incident on a longitudinally polarised 1 H 
gas target internal to the HERA storage ring at 
DESY. The beam polarisation was measured con- 
tinuously using Compton backscattering of circu- 
larly polarised laser light |l4|] . The average beam 
polarisation for the analysed data was 0.55. 

The HERMES polarised target Q consists 
of polarised atomic X H gas confined in a storage 
cell, which is a 40 cm long open-ended thin- walled 
elliptical tube located on the beam axis inside 
the HERA vacuum pipe. It is fed by an atomic- 
beam source of nuclear-polarised hydrogen based 
on Stern-Gerlach separation fl6|| . It provides an 
areal target density of about 7 x 10 13 atoms/cm 2 . 
The nuclear polarisation of atoms and the atomic 
fraction are continuously measured with a Breit- 
Rabi polarimeter and a target gas analyser |p^7|| . 
The average value of the target polarisation for 
the analysed data was 0.88 |18|. The fractional 
systematic uncertainties of the beam and target 
polarisations were 3.4% and 4.5%, respectively. 
The integrated luminosity for this data set was 
70 pb" 1 . 

Scattered positrons were detected by the HER- 
MES forward spectrometer, which is described in 
detail elsewhere Q . The kinematic requirements 
on the scattered positrons for the analysis in the 
nucleon-resonance region were: 1 GeV 2 < W 2 < 



applying data quality criteria, about 0.13 million 
events were selected. 

For all detected positrons the angular resolu- 
tion was better than 0.6 mrad, the momentum 
resolution was better than 1.6% aside from brems- 
strahlung tails and the Q 2 -resolution was better 
than 2.2%. The limited Fy 2 -resolution (about 840 
MeV 2 , or AW ~ 240 MeV) in the resonance re- 
gion did not allow the contributions of the in- 
dividual nucleon resonances to be distinguished. 
To evaluate the smearing corrections and the con- 
taminations intruding into the resonance region 
from the elastic and deep-inelastic regions, events 
were simulated using a Monte-Carlo code that in- 
cludes elastic, deep-inelastic and resonance contri- 
butions. The description of the resonance contri- 



bution was based on the model of Ref. 20 1 . The 
deep-inelastic region was modelled using the pa- 
rameterisation of Ref. |2l[] while the elastic form 
factors were taken from Ref. |§ . In Fig. | the dis- 
tribution of events as a function of W 2 is presented 
in comparison with the simulation. It is apparent 
that the shape of the simulated distribution agrees 
well with that of the data. It was found that the 
relative contaminations from the clastic and DIS 
regions in the yield of the resonance region range 
from 10% to 2% and from 7% to 16% respectively, 
as Q 2 increases from 1.2 GeV 2 up to 12 GeV 2 . 

Data were divided into six bins in Q 2 , but only 
one bin in W 2 . In each Q 2 -bin the average longi- 
tudinal asymmetry Ai for virtual photoabsorption 
was calculated using the formula 



r]A 2 



(6) 



where D and r\ are factors |13| that depend on 
kinematic variables. The quantity D depends also 
on R = olIot-, the ratio of the absorption cross 
sections for longitudinal and transverse virtual pho- 
tons. Ai is related to longitudinal-transverse in- 
terference. The cross section asymmetry A\^ is 
given by : 



N-L- 



N-Lp+N-Lp 



(7) 



4.2 GeV , 1.2 GeV < 



< 12 GeV . After 



Here, N^(N^) is the number of scattered posit- 
rons for target spin parallel (anti-parallel) to the 
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Figure 1: Comparison of the event distribution 
for W 2 > 1 GeV 2 (squares) with the Monte-Carlo 
simulation (histogram). An overall normalisation 
factor was applied to the simulation to match the 
data. Also shown are the smeared distributions 
from the elastic, resonance and the deep inelastic 
regions. The vertical lines indicate the resonance 
region considered in the analysis. 

beam spin orientation. The deadtime-corrected 
luminosities for each target spin state are L^*~' 

and Lp^' , the latter being weighted by the prod- 
uct of beam and target polarisations. The cross 
section asymmetry A|j es in the resonance region 
was corrected for contaminations originating from 
elastic and deep-inelastic scattering, as discussed 
above. The asymmetry for the elastic contribu- 
tion was taken from |22||, w hile for the DIS re- 
gion a parameterisation pa] based on world data 
has been used. Model dependent uncertainties due 
to these asymmetries and contributions from the 
Monte Carlo simulation are negligible, as this cor- 
rection is less than 5% for Q 2 < 5 GeV 2 . Radia- 
tive corrections were calculated using the codes 
described in Ref. [M, and were found not to ex- 
ceed 2% of the asymmetry AJ, es . The values for 



A\ es + r]A 2 es for the measured range of Q 2 values 
are presented in Table 1. The asymmetry A\ was 
evaluated using Eq. ^ under the assumption that 
A 2 = 0.06 |25| in the whole resonance region, and 
with an average depolarisation factor D weighted 
by the event distribution. 

The contribution I^fl) H of the resonance region 
to the GDH integral was determined in each Q 2 - 
bin from the asymmetry A\, according to Eq. |], 
in which the integration limits were determined by 
the 1 GeV 2 < W 2 < 4.2 GeV 2 range. The unpo- 
larised structure function F\ = F 2 (l-|-7 2 )/(2x(l + 
R)) was calculated from a modification of the pa- 
rameterisation of F 2 given in Ref. p(J that ac- 
counts for nucleon-resonance excitation assuming 
R = cxl/ut is constant and equal to 0.18 in the 
whole resonance region. It is worth noting that 
due to cancellation between the R-dependences 
of F\ and D at low y, the final result is insen- 
sitive to the choice of R. In the integration the 
W 2 dependence of the integrand F\jK within the 
individual bins was fully accounted for. 

The results for Tq^h are presented in Fig. 0. 
The integral strongly decreases with Q 2 over the 
entire measured range. The magnitude of the sys- 
tematic uncertainty is indicated by the band. The 
dominant contribution to the systematic uncer- 
tainty is due to the lack of knowledge of A 2 . This 
contribution (up to 15%) was evaluated from the 
total error quoted for a measurement in the reso- 
nance region: A 2 — 0.06 ± 0.16 [^5|. This range 
is consistent with two other possible assumptions 
for A 2 : A 2 = 0, or A 2 = 0.53MXA/Q 2 , which 
describes the behaviour in the deep-inelastic re- 
gion |ll|. Other contributions are uncertainties 
from the beam and target polarisations (5.3%), 
from the spectrometer geometry (2.5%), from the 
combined smearing and radiative effects (up to 
10%) and from the knowledge of F 2 (2%). The 
smearing contribution to the systematic uncertainty 
was evaluated by comparing simulated results from 
two very different assumptions for A\\ a power 
law (A\ — x - 727 ) that smoothly extends the DIS 
behaviour for the asymmetry into the resonance 
region Q, and a step function (A± = —0.5 for 
W 2 < 1.8 GeV 2 and A x = +1.0 for 1.8 GeV 2 < 
W 2 < 4.2 GeV 2 ) that is suggested by the hypothe- 
sis of the possible dominance of the P33-resonance 
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Figure 2: The GDH integral as a function of Q 2 for 
the region 1.0 GeV 2 < W 2 < 4.2 GeV 2 . The error 
bars show the statistical uncertainties. The white 
and the hatched bands represent the systematic 
uncertainties with and without the A 2 uncertainty 
contribution. The dashed M and the solid jl0| 
curves are predictions based on a Q 2 -evolution of 
nucleon-resonance amplitudes. 

at low W 2 and of the Su -resonance at higher W 2 . 

The results for I G e f) H are compared in Fig. || 
with two predictions for the contribution of nuc- 
leon-resonance excitation to the integral defined 
in Eq. |^. Burkert and Li j^] parameterised the ex- 
perimental Q 2 -evolution of the main nucleon reso- 
nances (P 33 (1232), P n (1440), Sn(1535), £>i 3 (1520), 
^15(1680)), and assumed single-quark transitions 
to evaluate the contributions from other resonances. 
Aznauryan jl0| described the resonance excita- 
tion in the approximation of infinitely narrow res- 
onances, and included a contribution from one- 
pion exchange in the near-threshold region. Both 
models predict a sudden decline in I t qdh as Q 2 
falls below 1.5 GeV 2 , due to a large negative con- 
tribution at low Q 2 by the helicity structure of the 
-P33-resonance. At higher Q 2 the P33-resonance 



magnetic form factor strongly decreases with in- 
creasing Q 2 , and the positive contribution to Iqjjh 
arising from the excitation of higher-mass reso- 
nances becomes dominant. Neither of these mod- 
els includes the non-resonant multi-hadron pro- 
duction channels, which should provide an addi- 
tional positive contribution for the region W 2 < 
4.2 GeV 2 . Comparison with the data suggests 
that for Q 2 ~ 1.5 GeV 2 , the resonance-excitation 
models are not sufficient to fully explain the exper- 
imental result for I G % H - Other predictions exist 
for the resonance-excitation contribution to gen- 
eralised GDH integrals, but they are limited to 
regions of lower Q 2 

To complete the evaluation of the full integral 
Igdh, data from the DIS region (4.2 GeV 2 < 
W 2 < 45 GeV 2 ) were reanalysed in the same Q 2 - 
bins as for the kinematically more restricted reso- 
nance region, following the procedure described in 
a previous HERMES publication Q. A total of 
1.52 million events were selected in this IV 2 -range. 
The systematic uncertainty for this region is the 
same as published in |ll|[l^]. The systematic un- 
certainty on A2 in DIS region does not contribute 
significantly. 

In Table 1 the resonance region contribution 
IqI> h , the integrals in the full measured 

region and the full GDH integrals are reported. 
The latter was calculated in each Q 2 -interval by 
adding to an estimate of the unmeasured 

contribution for W 2 > 45 GeV 2 . This was calcu- 
lated using a multiple-Reggeon exchange param- 
eterisation || for ai(v,Q 2 ) — a3.(v,Q 2 ) at high 
energy, and amounted to about 3.5 fib for all Q 2 - 
bins. A parameterisation for g\ [^J based on a 
NLO-QCD analysis provided within 5% the same 
results as the multiple-Reggeon exchange analysis. 
This difference was taken as the systematic uncer- 
tainty of the high energy contribution. It is worth 
noting that for the low Q 2 -bins, both the statisti- 
cal and the systematic uncertainties of Igdh are 
dominated by those from the resonance region. In 
this region, these uncertainties are large due to 
the smallness of D and to the large size of 77, re- 
spectively. 

In Fig. ||a) the total GDH integral is shown to- 
gether with the partial integrals for W 2 < 4.2 GeV 2 
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Figure 3: a) Igdh as a function of Q 2 for various 
upper limits of integration: W 2 < 4.2 GeV 2 (tri- 
angles), W 2 < 45 GeV 2 (squares), and the total 
integral Igdh (circles). The squares have been 
slightly shifted to make them more visible. The 
curve is the Soffer-Teryaev model || for the total 
integral, b) IgdhQ 2 j (167r 2 a) as a function of Q 2 . 
For both panels, the error bars show the statisti- 
cal uncertainties, and the white and the hatched 
bands at the bottom represent the systematic un- 
certainties for the total integral with and without 
the A2 contribution. 

and for W 2 < 45 GeV 2 . The contribution of the 
resonance region to the full GDH integral is small 
for Q 2 values above 3 GeV 2 . 

Fig. |3|a) also shows a prediction (8j based on a 



Q 2 -evolution of the two polarised structure func- 
tions gi and g2, without consideration of any ex- 
plicit nucleon-resonance contribution. This pre- 
diction is in good agreement with the experimen- 
tal data. 

In the whole energy range, Igdh is consistent 
within the uncertainties (x 2 /^df = 0.4) with a 
simple jJQ 2 power law. This is demonstrated 
in Fig. ||b where the results for Igdh are multi- 
plied by Q 2 /(167r 2 a). In the deep-inelastic limit, 
this quantity is equivalent to Ti (see Eq. ||). The 
present results are in agreement with the measure- 
ments of Ti = 0.120 ± 0.016 at Q 2 = 10 GeV 2 @ 
and Ti = 0.129 ± 0.010 at Q 2 = 5 GeV 2 [||. In 
addition, values of Ti extracted from the present 
data are also consistent with a measurement of 
Ti = 0.104±0.017 at Q 2 = 1.2 GeV 2 |§ in which 
the structure function g\ was measured in the res- 
onance region. 

In summary, the Q 2 -dependence of the gener- 
alised Gerasimov-Drell-Hearn integral for the pro- 
ton was determined for the first time in both the 
resonance and the deep-inelastic regions, covering 
the Q 2 -range from 1.2 to 12 GeV 2 . In the res- 
onance region, the data suggest that for Q 2 > 
1.5 GeV 2 , existing resonance-excitation models are 
not sufficient to fully explain the experimental re- 
sult for Iffi H . Above Q 2 = 3 GeV 2 the DIS con- 
tribution to the generalised GDH integral is dom- 
inant. The Q 2 -behaviour of Igdh suggests that 
there are no large effects from either resonances 
or non- leading- twist, and indicates that the sign 
change of Igdh to meet the real photon limit oc- 
curs at Q 2 lower than 1.2 GeV 2 . 

We thank S. Gerasimov, V. Burkert and I.G. 
Aznauryan for useful discussions and I.G. A. for 
providing the curves of her calculations. We grate- 
fully acknowledge the DESY management for its 
support, the staffs at DESY and the collaborat- 
ing institutions for their significant effort, and our 
funding agencies for financial support. 
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s 



Table 1: Results for A\ es + r]A T 2 es , the resonance part(/£fg H ) to the GDH integral, and the total measured 
integral {Iq^h), as well as the full GDH integral (Igdh), including the unmeasured part . Errors represent 
the statistical uncertainty for A\ es + rjA r 2 es , and the statistical and the systematic uncertainties of the 
integrals. 



< Q 2 > 


A\ es + r)A r 2 es 


rres 
GDH 


Tineas 
GDH 


Igdh 


[GeV 2 ] 




[ M b] 


[Mb] 


[ M b] 


1.5 


0.71±0.16 


21.4±5.2±4.1 


37.9±5.5±5.1 


41.2±5.5±5.1 


2.1 


0.77±0.18 


10.3±2.5±2.0 


24.3±2.8±2.9 


27.8i2.8i2.9 


2.7 


0.74±0.22 


4.9±1.5±0.9 


17.5±1.8±1.8 


21.0il.8il.8 


3.5 


0.79±0.22 


2.4±0.7±0.4 


13. Oil. Oil. 2 


16.5il.0il.2 


4.5 


0.97±0.29 


1.3±0.4±0.2 


8.8i0.6i0.8 


12.3i0.6i0.8 


6.6 


0.55±0.23 


0.08±0.03±0.01 


5.3i0.3i0.5 


8.6i0.3i0.5 
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